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NUMERICAL  COMPUTATION  OF  INVISCID  FLOW  AROUND  NOSE  TIP  OF 
SPACE  SHUTTLE  WITH  NND  SCHEME 

ABSTRACT:  The  NND  scheme  presented  by  Professor 
Zhang  Hanxin  of  the  CARDC  has  been  used  in  this  paper  to 
complete  the  following  computations;  the  two-dimensional 
inviscid  flow  around  the  symmetric  plane  of  the  nose  tip 
of  the  space  shuttle;  the  three-dimensional  inviscid  flow 
around  the  nose  tip  of  the  space  shuttle;  the  three 
dimensional  inviscid  flow  around  the  axisymmetric  indented 
body.  Under  the  general  coordinate  system,  the 
computation  grids  have  been  generated  by  algebra 
equations.  The  characteristic  compatibility  equations 
have  been  used  to  treat  the  shock  and  body  surface 
boundaries,  and  there  is  a  special  treatment  in 
translation  of  the  coordinate  system.  Finally,  the 
comparative  computations  have  been  completed  using  second 
order  SCM  schemes.  The  results  of  computation  illustrate 
that  the  NND  scheme  used  in  this  paper  can  avoid  numerical 
oscillations  and  high  resolution  to  shock-capturing. 


Since  J.V.  Rakich^'^  used  MacCormack  two  step  grid  computation 
of  the  space  shuttle  flow  field,  several  people  have  done  work 
in  this  area.  For  the  mixed  zones  of  the  nose  cone,  the 
stabilization  method  was  used.  For  the  super  sonic  areas  of  the 
body  of  the  space  craft,  the  propulsion  method  was  used.  For  the 
external  shock  wave,  an  assembly  method  was  used.  These  methods 
are  the  general  methods  for  such  complex  shapes  as  a  space  shuttle. 
However,  engineers  have  always  been  searching  for  a  simple 
calculation  method  with  a  small  amount  of  computation  which  is 
deeply  rooted  in  physics. 

In  1983,  A.  Harten  proposed  a  total  variation  decrease 
scheme®^.  This  scheme  is  very  effective  in  capturing  the  shock 
wave  and  can  also  eliminate  the  normal  two  step  scheme  for 
capturing  the  vibration  generated  when  capturing  the  shock  wave. 
Later,  this  scheme  based  on  the  idea  of  flux  limitation  was  widely 
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expanded^"^®^ .  In  1985,  H.C.  Yee  used  the  TVD  scheme  to  do 
computations  for  the  NACA  0012  wing  form  on  a  graphic  coordinate 
system,  and  S.R.  Chakravarthy  applied  the  TVD  scheme  to 
computations  for  complex  shapes,  both  with  good  results. 

Last  year,  at  the  discussion  section  on  "The  TVD  Scheme  and 
Its  Applications"  held  by  the  graduate  department  at  the 
Aerodynamics  Research  Center,  Professor  Zhang  Hanxin  proposed  a  NND 
scheme  which,  when  compared  to  other  TVD  schemes,  has  the  advantage 
of  being  simple  and  having  a  small  amount  of  computations.  This 
article  will  use  this  NND  scheme  to  compute  the  two  dimensional 
inviscid  flow  around  the  symmetrical  cross  section  of  the  nose  tip 
of  the  space  shuttle  and  the  three  dimensional  inviscid  flow  around 
the  nose  tip  of  the  space  shuttle  and  the  axisymmetric  indented 
body.  The  results  indicate  that  it  is  satisfactory  for  both  the 
scheme  convergence  speed  and  the  capture  of  the  shock  wave  within 
the  flow  field. 
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II:  CONTROL  FORMULA  AND  BOUNDARY  CONDITIONS 

1 .  CONTROL  FORMULA 


For  a  rectangular  coordinate  system, 
written: 


Euler’s  formula  may  be 
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We  obtain  the  control  formula  for  the  object  surface  coordinates 


In  this  formula: 


.  AIL.  *'<7  -0 

On 


If: 


U-U/J 

-(I.  a  4  c.  E 

'  F  (r),  -i  77.  h  -  F  J 

G  -  (C,  u  s\  E  f  C,  F  r  C.  f?)  / ; 


OG 

'^U 


Then  Tk,  S,  C  have  the  characteristic  values  of: 

A,  — n,—  n.ti  r  f!,v —  C\/  n\  +n}  +nl 
' ..~n,  +  n,  u-rn,v-rn,iu 

=  n,  -r  n,  u4-n,i;  +  «,iw  +  Cv/  n;  -r  nj  f-  nl 

n  represents  4  ,  jj  ,  or  ^  . 
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2 .  BOUNDARY  CONDITIONS 

On  the  surface  of  an  object,  (  rj  ^0  )  ,  a  condition  for 
impermeability : 

On  a  shock  wave  surface,  {  =  using  the  R-H  relationship 

in  a  movement  coordinate  system,  we  can  obtain: 
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Herein,  the  subscript  "2"  represents  the  post  wave  value  and  N  is 
the  vector  of  the  front  of  the  shock  wave. 

Ill:  DIFFERENCE  SCHEME  AND  BOUNDARY  PROCESSING 

1 .  DIFFERENCE  SCHEME 

This  article  uses  the  non  division  explicit  difference  scheme 

u::i -•> 

'  -.^G  . O'  <  ;) 
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Herein:  A;  ^! 

"•‘"'aT’  A,,  ’  Ay 

'.T-/:;  .1,). 
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A  similar  expression  can  also  be  written  for  f,.' 

0  i.y<0 

nunmod(.v,  y)—  '  x  |xi<|y. 

:j  \x\>\y 

which  are  called  limiting  factors. 

2 .  BOUNDARY  PROCESSING 

Because  the  boundary  conditions  of  the  physical  object  and  the 
boundary  conditions  of  the  shock  wave  are  incomplete,  it  is 
necessary  to  have  four  additional  formulas  for  the  surface  of  the 
object  and  one  additional  formula  for  the  shock  wave.  According  to 
characteristic  theory,  in  fluid  motion,  the  information  is 
broadcast  along  a  characteristic  direction.  On  a  shock  wave,  there 
is  a  characteristic  direction  orientated  with  the  shock  wave.  The 
formula  for  its  consistent  relationships  can  be  used  to  supplement 
shock  wave  conditions.  On  the  surface  of  the  object,  there  are 
four  characteristic  consistent  relationship  formulas  which  can  be 
used  to  supplement  object  surface  conditions.  In  this  manner,  the 
object  surface  and  the  shock  wave  boundary  formulas  are  confined. 
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On  the  shock  wave,  ;.,>0  •  Its  consistent  relationship  formula 
can  be  written  as: 
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After  we  obtain  the  post  wave  pressure  p2,  we  insert  it  into 
the  shock  wave  boundary  conditions,  and  we  an  obtain  other  post 
wave  parameters . 

Because  the  object  surface  is  static,  the  impermeable 
conditions  can  be  written  as: 


Ou 
'  dr 
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On  the  object  surface,  -'0  .  its  characteristic  consistent 
relationship  formula  is  substituted  for  by  the  object  surface 
impermeable  conditions.  After  insertion  into  the  impermeable 
conditions,  the  object  surface  value  boundary  conditions  are: 
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sr 

1:,  H.  S: 

s-  2'/.//.  -  s, 

5  H.  S. 

0  t> 

Subscript  ij  a  vector  matrix  element  of  the  B  left  characteristic 
matrix. 

By  collating  lambda,  we  obtain: 

=-v'2  C(R,-S,) 

<1^  '  J 2  ^  ^  •  7  •  S' )  •  n  -  s  ■: 

>-  -  S,)] 

4v  =  p  ^.(A-  +5.)] 

'  7,'/-.  -S',)]  7  ».,..  =  7..,.,  /  v/Vf-  7:  -  r;; 

From  the  formulas  above  we  can  obtain  the  object  surface  physics 

quantities . 

At  the  object  axial  boundary,  using  span  axial  series 
differential  equations  we  can  avoid  the  irregularities  produced  by 
three  dimensional  coordinate  conversion.  For  symmetric  surfaces, 
we  use  symraetr-y  to  carry  out  parameter  extrapolation  of  the 
supplementary  boundary  values. 

IV:  COMPUTATION  EXAMPLES 

1.  TWO  DIMENSIONAL 

In  order  to  check  the  stability  and  convergence  rate  of  the 
differential  scheme  as  well  as  the  suitability  of  the  boundary 
processing,  this  article  first  applies  the  scheme  to  space  shuttle 
nose  tip  symmetrical  cross  section  two  dimensional  non  viscous  flow 
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computations.  In  order  to  speed  up  the  convergence,  it  uses 
partial  time  steps. 


i'  \  !  (F;  -i  Cv/'  t)]  ■  ’ll 

In  this  formula,  U  and  V  are  velocity  along  the  directions  xi  and 
eta.  Delta  time  subscript  reference  is  determined  by  the  numerical 
test.  The  graph  was  61  by  11,  and  the  initial  shock  wave  form  was 
a  parabola. 


ILLUSTRATION  ONE;  ISOBARIC  CHART  (NND  SCHEME) 

Illustration  one  is  the  isobaric  chart  of  the  flow  field  when 
the  angle  of  attack  is  zero.  Illustration  two  is  a  comparison  of 
the  pressure  on  the  object  surface  obtained  through  computation 
under  conditions  of  viscous  and  non  viscous  flow.  Illustration 
three  is  a  comparison  of  the  computed  results  using  the  NND  and  the 
SCM  schemes.  From  these  illustrations  we  can  see  that  the  NND 
scheme  not  only  can  do  a  good  job  in  capturing  the  shock  wave  in 
the  flow  field,  but  does  not  generate  oscillation  near  the  shock 
wave.  From  Illustration  two  we  can  see  that  because  of  the 
existence  of  separated  areas,  the  shock  wave  in  the  viscous  flow 
calculations  is  much  weaker  that  in  the  non  viscous  flow. 
Therefore,  changes  in  object  surface  pressure  tend  to  be  more 
gradual . 
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ILLUSTRATION  TWO:  COMPARISON  OF  OBJECT  SURFACE  PRESSURE 
WITH  VISCOUS  AND  NON  VISCOUS  FLOW 

NOTE:  Solid  dots  are  viscous  flow,  and  hollow  dots  are  non 

viscous . 


ILLUSTRATION  THREE;  COMPARISON  OF 
SURFACE  PRESSURE  DISTRIBUTION 


2 .  THREE  DIMENSIONAL 


Three  dimensional  computations  include  two  parts, 
shuttle  nose  tip  and  the  axisymmetric  indented  body, 
number  is  10.1  and  no  consideration  is  given  to  yawing, 
step  is; 
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The  graph  is  41  y  11  by  17.  The  initial  shock  wave  form  is 
presented  as  a  spiral  parabolic  surface.  The  space  shuttle  nose 
tip  calculated  angle  of  attack  was  calculated  at  both  0.0  and  minus 
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5.0  degrees.  The  indented  body  calculations  were  made  for  both  0.0 
and  minus  5.0  degrees. 

Illustration  four  is  the  calculated  result  of  the  space 
shuttle  nose  tip  flow  field  at  an  angle  of  attack  of  minus  5.0 
degrees.  Illustration  five  is  the  meridian  surface  pressure 
distribution  calculated  with  the  NND  scheme  and  the  SCM  scheme 
compared  with  as  bibliographical  item  number  19.  The  results  are 
fairlv  consistent. 


..'f  5  10 ,  a  =  -  5  0  .NND.  ;M  =  lO ,  a  =  -  o.O  ,N‘ND  (M  -  1' .  u  =  -  i.C' .  NND 


ILLUSTRATION  ONE  ILLUSTRATION  TWO  ILLUSTRATION  THREE 
ISOBARS  ISOPYCNAL  ISOMACH 


ILLUSTRATION  5.1:  COMPARISON  OF  OBJECT  SURFACE 
PRESSURE  DISTRIBUTION 

NOTE;  Solid  dot  is  NND  scheme.  Hollow  dot  is  SMC  scheme. 
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ILLUSTRATION  5.2:  COMPARISON  OF  OBJECT  SURFACE 
PRESSURE  DISTRIBUTION 

NOTE:  Solid  dot  is  NND  scheme,  hollow  dot  is  SCM  scheme. 


ILLUSTRATION  SIX 
ISOBARIC  GRAPH 


Illustration  six  is  an  isobaric  graph  of  the  indented  body 
flow  field  at  an  angle  of  attack  of  five  degrees.  Illustration 
seven  is  pi  meridian  object  surface  pressure  distributions 
calculated  with  the  NND  and  SCM  schemes  compared  with  bibliographic 
item  19.  The  results  also  conform  very  well. 
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ILLUSTRATION  7.1:  OBJECT  SURFACE  PRESSURE  DISTRIBUTION  COMPARISON 


ILLUSTRATION  7.2;  OBJECT  SURFACE  PRESSURE  DISTRIBUTION  COMPARISON 

V:  CONCLUSIONS 

1.  The  NND  scheme  used  in  this  article  is  not  only 
structurally  simple  and  has  a  small  amount  of  computations,  its 
results  in  computing  complex  forms  are  also  satisfactory. 

2.  The  NND  scheme  used  in  this  article  can  do  a  good  job  of 
capturing  shock  waves  in  flow  fields,  it  has  fairly  high 
discrimination,  and  it  does  not  create  physical  vibrations  near  the 
shock  wave. 

3.  In  the  production  of  a  three  dimensional  graph,  this 
article  used  span  axial  series  differential  equations,  thus 
avoiding  irregularities  produced  by  coordinate  conversion,  and 
reducing  the  amount  of  computation.  The  computational  graph 
generated  using  algebra  formulas  is  not  only  quick  with  a  small 
amount  of  computation,  but  it  also  is  accurate  and  meets 
requirements . 
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